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The one-dimensional Kondo lattice model with attractive interaction among the conduction elec- 
trons is analyzed in the case of half-filling. It is shown that there are three distinct phases depending 
on the coupling constants of the model. Two phases have a spin and charge gap. While one shows 
a clear separation of the spin and charge excitation spectrum the other phase may be characterized 
as a band insulator type where both excitations are due to two-particle states. The third phase 
is gapless in both channels and has quasi long-range order in the spin and charge density wave 
correlation. In this phase the spin and charge excitations have again a clearly separated spectrum. 
For the analysis we discuss first two limiting cases. Then a density matrix renormalization group 
calculation on finite systems is applied to determine the phase diagram and the correlation functions 
in the gapped and gapless phase for general couplding constants. 
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I. INTRODUCTION 

The Kondo lattice model (KLM) describes a many- 
body system of two distinct types of degrees of freedom, 
itinerant electrons and localized spins which are arranged 
on a regular lattice. This model can be considered as an 
extension of the single impurity Kondo model where elec- 
trons interact with a single localized spin. The Hamilto- 
nian of the KLM consists of two parts, the kinetic energy 
of the itinerant conduction electrons and the local ex- 
change interaction between electron spin Sc and localized 
spin Sf, both spin 1/2 degrees of freedom. 



(1) 



where the operator Cis (c|^) annihilates (creates) a con- 
duction electron on site i with spin s (=Tji) {S^^ — 
{fi/2) J2s s' '^is'^ss''^is') and the sum in the first term runs 
over all nearest neighbor bonds Furthermore, t 

denotes the hopping matrix element and J is the antifer- 
romagnetic exchange couplingtJ. 

In this model both the conduction electrons and the 
localized spins are each uncorrelated. Correlation only 
appears through the exchange interaction. It is impor- 
tant to realize that the understanding of the single im- 
purity Kondo problem does not simply extend to the lat- 
tice case. Indeed the research during recent years has 
shown that complicated correlation effects occur within 
the lattice model yielding a variety of physical phenom- 
ena beyond the single impurity picture. Although the 
KLM is certainly of interest of its own as a generic 
model of strongly correlated electrons, studies of this 
model are also motivated by real materials such as the so- 
called heavy fermion compounds or the Kondo insulators. 
Both systems where correlation effects between itinerant 
and localized electrons clearly dominate the low-energy 
physics. 



In this article we discuss the extension of the KLM 
obtained by including the direct interaction among the 
conduction electrons. We use the most simple form for 
the interaction, which acts only between electrons on the 
same site, 



Hi, 



U 
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We call Ti. = T^klm + ^int Kondo-Hubbard model 
(KHM) . This model was recently also considered in con- 
nection with the ferromagnetic ground state away from 
half fillinga. In the following we will focus on the one- 
dimensional (ID) system with a half-filled electron band. 

The half-filled KLM is considered as a good starting 
point to understand the Kondo insulators, a class of ma- 
terials which show a spin and a charge gap at low tem- 
peratures. In contrast to ordinary band insulators the 
two gaps are different, indicating a separation of the spin 
and charge degrees of freedom due to correlation effects. 
The half-filled KLM shows indeed this type of proper- 
ties. The properties of the spins are dominated by short- 
ranged antiferromagnetic correlations, i.e. this state can 
be considered as a spin liquid. A particular feature of 
the ID KLM is that the spin liquid state exists for all 
finite values of J. This could be shown numerically us- 
ing exact diagonalizationa or, more recently, the density 
matrix renormalization group tedrmiquoj and the map- 
ping to a non-linear sigma modelEI. (Note that in higher 
dimensions a transition between the spin liquid and an 
antiferromagnetically ordered state is expected at a crit- 
ical value of J. This has, however, not been established 
so far either by analytical or numerical methods.) For 
weak coupling, J <C i, one finds that the spin gap (Ag) 
and charge gap (Ac) depend in a very different way on 
J, 



As oc exp(— 1/apJ) 
Ac cx J 



(3) 
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where p is the density of states at the Fermi level for the 
free electrons. The spin gap energy gives an energy scale 
formally related to the Kondo temperature Tk of a single 
localized spin (a — 1). For the lattice of localized spins 
a is enhanced by aj/alue of about 1.4 obtained from nu- 
merical simulationsO. The similarity with the Kondo en- 
ergy scale indicates that the ground state of the half-filled 
KLM corresponds to a singlet bound state like the Kondo 
singlet state. However, correlation effects among the lo- 
calized spins tend to increase the binding energy through 
the formation of a collective singlet state. On the other 
hand, the charge gap has no counter part in the single im- 
purity case (where the system forms a Fermi liquid). The 
charge gap proportional to J originates from the strong 
antiferromagnetic correlation of the localized spins. Al- 
though short ranged it provides a staggered background 
for the electron motion, which yields the features ofj-a 
doubled unit cell as found in a spin density wave stateQ. 
Thus the energy scale is set by the exchange coupling be- 
tween the electrons and the localized spins. Also in the 
strong coupling limit (J ^ t) we find a clear distinction 
between the two excitations as we will consider below. 

The spin liquid state of the KLM can be characterized 
as the formation of a collective Kondo singlet involving 
the conduction electrons and localized spins. The sin- 
glet formation is optimized if configurations of conduc- 
tion electrons with doubly occupied and empty sites are 
suppressed. Both configurations "remove" the electron 
spin degree of freedom so that the exchange term can- 
not be active. Therefore strong charge fluctuations of 
the conduction electrons tend to weaken the spin liquid 
state. The inclusion of a repulsive interaction between 
the conduction electrons as in Eq.(2) with [/ > leads 
to a suppression of the charge fluctuations. It is a well- 
known fact that the Hubbard model with repulsive inter- 
action develops a charge gap at half-filling for any finite 
U in one dimension, while the spin excitation remain gap- 
less. Consequently, we expect that the spin liquid state 
is further stabilized by the repulsive interaction. This is 
clearly seen in numerical calculations where, in particu- 
lar, in the limit oi J <^t the spin gap is enhancecLthrough 
an increase of the factor a in the exponent Eq.(3)El. While 
the spin gap goes to zero for J ^ the charge gap re- 
mains finite, if C/ > 0. 

We now ask what will happen if the interaction among 
the conduction electrons is attractive. There is a relation 
between the positive and negative U Hubbard model due 
to particle-hole symmetry at half filling. The change of 
sign for U leads effectively to an exchange of the charge 
and spin degrees of freedom. Indeed the charge excita- 
tions can be described as isospins^completely analogous 
to the spins, as we will see belowQ. For U < the spin 
sector of the conduction electrons has a gap while the 
charge (isospin) excitations are gapless. Obviously this 
spin gap weakens the spin liquid phase characterized by 
the formation of the singlets between localized and elec- 
tron spins. Therefore a competition arises between the 
the spin and the charge fluctuations whose behavior is 



determined by the relative strength of the coupling con- 
stants J and U . 

It is the goal of this paper to investigate this compe- 
tition for the case of C/ < 0. By the analysis of two 
limiting cases we will demonstrate that the attractive in- 
teraction leads to a phase transition where the character 
of the ground state and the excitations change qualita- 
tively. The state which is in competition with the spin 
liquid phase may be characterized by the property that it 
has quasi-long-range order in the spin and charge sectors. 
The dominant correlations are that of an antiferromag- 
net for the localized spins and of a charge density wave 
for the conduction electrons, i.e. a spin-charge density 
wave (SCDW). Additionally we find that within the spin 
liquid a phase can occur where the spin and charge ex- 
citations have the same energy scale and behave similar 
to a band insulator. In the following we will analyze 
these properties first by characterizing the states in the 
limit t ^ |[/|, J. In a second step we will investigate the 
phase transition and the phase diagram, U versus J, by 
means of numerical simulation using the density matrix 
renormalization group method (DMRG). 

II. PHASES IN TWO LIMITING CASES 

In this section we show that there are two distinct 
phases for the KHM with attractive interaction. To this 
end it is helpful to consider two limiting cases which al- 
low a simple analysis of the model. These limits are 
t < J < ~U and t < -U < J. At the very begin- 
ning let us start with t = 0, the atomic limit. Then 
the states can be represented most easily in a real space 
basis. For U — the ground state j^E's) is a product 
of onsite singlets, i.e. on each site we find one conduc- 
tion electron which forms a spin singlet with the localized 
spin (l^fs) = {cl^\ ii) - c^J ti))/V2; \si) is the electron 
vacuum with the localised spin Si on the site i). This 
state is non-degeneratm. The lowest spin excitation |5't) 
{— cl^\ Ii)) corresponds to the spin triplet configuration 
on one site. Since one singlet had to be replaced by a 
triplet the excitation energy is J. This state is highly 
degenerate (3A^) due to the freedom of position and the 
three triplet orientations {N : number of lattice sites). 
The lowest charge excitation j^fc) consists of a doubly- 
occupied and an empty site, which we call doublon and 
holon, respectively (|5'c) = + c'-^c'-^\sj)). The excita- 
tion energy is 3 J/2, because two singlets are destroyed to 
generate this state. The doublon and holon are fermionic 
particles with a spin 1/2 degree of freedom, which com- 
bine to a spin singlet for a pure charge excitation. Then 
the degeneracy is N{N — 1). The combination of spin 
and charge configuration corresponds to the triplet con- 
figuration of the doublon-holon spins with a degeneracy 
of 3N{N -1). 
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FIG. 1. Level scheme of the states for t — as & function 
of the electron-electron interaction U. 

Turning on U now we can change the level scheme 
(Fig. 0). The relative position of the singlet and triplet 
states, I^E's) and |^t), is unchanged. However, the 
charge (doublon-holon) excitation is shifted according to 
3 J/2 + U. For positive U no qualitative change occurs 
in the level scheme. Negative U, however, lead to a re- 
arrangement of the onsite energy levels for sufficiently 
large \ U\. For U = — J/2 the charge excitation passes the 
triplet state and for U — —3 J/2 the singlet state. This in- 
dicates that the attractive interaction between electrons 
yields a qualitative modification of the system. Now let 
us discuss the situation where t is finite and the above 
mentioned degeneracies are lifted. 



A. The spin liquid state: \U\ <C J 

For small t we can now use perturbation theory to 
describe the effect of hopping of the conduction electrons. 
A detailed discussion of this type of perturbation can be 
found in Ref. ||. For the limit Li/| ^ J the singlet state 
l^'s) remains the ground stateu. Its energy acquires a 
correction in second order due to the polarization of the 
onsite singlet, i.e. the neighbor electron and localized 
spin are involved in the formation of the singlet which 
becomes more extended. 

The hopping term leads to the mobility of the spin 
triplet and the doublon-holon excitations. The former 
behaves as a single quasiparticle with an excitation en- 
ergy 



Etiq) = J 



At' 



3J + 2U J + 2U 



(1 — cos 5) 
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where q is the momentum (lattice constant a = 1). The 
excitation energy for the doublon-holon state has the 
form of two-particle excitation with an energy depend- 
ing on two momenta k and q, 



3J 

T 



u 



3J + 2U 
+t{cos{k + q) + cos(fc)} 



J 



(5) 



3J + 2U 



{cos(2(fc + (7)) +cos(2fc)}. 



Thus, the doublon-holon excitations form a continuum. 
The spin triplet excitation may be understood also as a 
bound state of a doublon and holon with their spins in the 
triplet configuration. In this sense the triplet excitation 
has the character of an "exciton" within the gap between 
the singlet ground state and doublon-holon continuum. 

The discussion of the level scheme for t = suggests 
that there is some change in the properties of the spin 
liquid phase as we turn U from towards negative values. 
There is a critical value U = — J/2 where the doublon- 
holon state falls below the triplet excitation (for t — 0). 
If U is smaller than — J/2 the triplet excitation discussed 
above is absorbed into the continuum of the spin triplet 
channel of the doublon-holon continuum. Therefore the 
lowest spin and charge excitations have both two-particle 
character like a particle-hole excitation. This state has 
still a finite gap and the structure of the excitations is 
essentially similar to that of a band insulator. The quasi- 
particles, the doublons and holons, are spin 1/2 fermions, 
composed of conduction electrons and localized spins. 

At the second critical value of U (= —3 J/2) the holon- 
doublon state passes the singlet ground state. Here the 
character of the ground state has to change so that this 
critical value should belong to a phase transition. In the 
following we want to discuss the properties of the new 
state. 

The perturbative treatment of the different states 
given above does not allow the evaluation of the tran- 
sition points for finite t. Close to the transition points 
a more complicated perturbation for degenerate states 
would be necessary. This leads, however, already to a 
complicated many-body problem. Thus, we leave the 
discussion of the phase boundary lines to the numerical 
part of our paper. 



B. The SCDW state: \U\ > J 

In the limit |C/| 3> J the ground state is highly degen- 
erate if we assume t — 0. A complete basis of these de- 
generate states is given by all real-space configurations of 
holons and doublons, with an equal number of each (half 
filling). Both onsite singlets jvPg) or triplets |\l/t) are now 
excited states which always should occur in pairs, due 
to the conservation of the electron number, with energies 
Ess = -3J/2-J7, Est - J/2-U andEtt J/2-C/for 
a singlet-singlet, singlet-triplet and a triplet-triplet pair. 
All these states are highly degenerate. 
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In order to lift the degeneracy of the ground state we 
include now the hopping process, assuming a small, but 
finite t. We can now generate an effective Hamiltonian in 
the Hilbert-subspace of doublons and holons only. At this 
point it is convenient to introduce the notion of isospin 
I (1=1/2) for the charge degrees of freedom of the holon 
and doublon. We identify a holon with isospin up and a 
doublon with isospin down. Like the ordinary spin also 
the isospin transforms according to the SU(2)-symmetry 
group. Additionally we introduce a phase convention by 
dividing the lattice into two sublattices A and B. Then 
the real space basis functions should be multiplied by 
n,e_Bexp(i7r4). 



(6) 



where Si = +1 ii i £ A and Si — —1 ii i £ B. With 
this notation the effective Hamiltonian which lifts the 
degeneracy in lowest order perturbation theory has the 
form, 



Weft = ^ 



I, I, 



[aS/j • Sfj + b] . 



(7) 



where the two constants are 
2^2 t 



u-i 



b = - 



MV4 3*2/2 tV4 



(8) 



u-i 



U + i 



u + M- 



(derived in the Appendix). In this formulation the S0(4)- 
symmetry of both the spin and isospin is obvious (SO (4) 
= SU(2) X SU(2)). The spin degrees of freedom are only 
due to the localized spins while the conduction electron 
spin is completely quenched in the effective model. Thus 
the conduction electrons only appear as charge degrees of 
freedom, i.e. isospins. It is easy to see that the coupling 
constant a is negative for U < — 3J/2 and vanishes for 
J ^ 0, because this leads to a complete decoupling of the 
localized spins from the conduction electrons. Note that 
the positive constant b determining the spectrum of the 
conduction electrons remains finite for J = so that the 
isospin degrees of freedom remain coupled. The electron 
spins appear only in states including onsite singlet or 
triplet components which are much higher in energy. 

For this type of Hamiltonian the following two exact 
statements are known. (1) For a finite system we can 
show that the ground state is a singlet in both the spin 
and the isospin channel and it is non-degenerate. This 
can be proven using a generalization of Marshall's theo- 
rem. (2) In the thermodynamic limit the excitations in 



both channels are either gapless or the ground state is de- 
generate with spontaneously broken parity. This can be 
demonstrated using a generalization of the LigttSchultz- 
Mattis theorem as given by Affleck and LiebElEl. There 
is no need to reproduce the proofs of these two state- 
ments here, since they are completely analogous to those 
applied to the Heisenberg spin system. 

The latter statement proves that in the limit 3> J 
the phase is different from the spin liquid state. While 
the spin liquid phase has a unique ground state and a 
gap in all excitations, this phase has either a degener- 
ate ground state or gapless excitations. In the case of 
gapless excitations we expect that the system has quasi- 
long-range order in both spin and isospin analogous to 
the ID Heisenberg model. The dominant correlation of 
both degrees of freedom is "antiferromagnetic" , which 
corresponds to the usual charge density wave correlation 
for the isospin part. Due to the qualitative difference of 
the two ground states, we conclude that there must be a 
phase transition between the two limits, depending on U 
and J. In the following we will use numerical methods 
to demonstrate that the state for U <^ — J has actually 
gapless excitations. 



III. THE PHASE DIAGRAM BY DMRG 

In order to determine the transition line in the phase 
diagram, U versus J, and the character of the different 
phases, we calculate the excitation energies and the cor- 
relation functions of the ground state in the KHM nu- 
merically. For this purpose it is necessary to treat long 
enough systems such that the transition between differ- 
ent phases can be observed reliably. We use the density 
matrix renormalization group (DMRG) method which al- 
lows us to study long chains by iteratively enlarging the 
system size and obtain ground, state wave functions with 
only small systematic errors. Ilil 

We calculate the elementary excitations for various U 
to identify the three phases discussed in Sec. II. The 
DMRG scheme is designed to obtain a very good approxi- 
mation of the ground state of a model. For the evaluation 
of excited states and their energies we consider the model 
in those Hilbert-subspaces which contain these states as 
lowest energy states. The real ground state is found for 
the subspace with S^^^^ = and = L {Nci number of 
conduction electrons). The spin excited state is obtained 
for S'^Qj — 1 and Nc = L and the charge excited state 
for S^^^ — and — L ± 2. In these sectors the low- 
est energy state is calculated for various system sizes up 
to L = 48 with the finite system algorithm using open 
boundary conditions. The extrapolation to the bulk limit 
is obtained from the scaling laws. For the gapped spin 
liquid phase the form A(L) = A(oo) + /3L~2 + 0{L-^) 
is assumed, since the lowest excited state generally cor- 
responds to the bottom of an excitation spectrum which 
can be expanded in terms of the square of the momentum. 
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fc^. On the other hand, for the region which is supposed 
to possess the gapless phase the bulk hmit value of the 
excitation energies is estimated simply using 1 /L scaling 
and 1 / log L scaling. 
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FIG. 2. Excitation energies of the half-filled Kondo lattice 
model with attractive interactions. Typical truncation errors 
in the DMRG calculations are 1 x 10"*^ for U/J = -1.4. Gap 
energies, exchange constant J, and Coulomb interaction U 
are in units of t. 
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FIG. 3. Excitation energies of the half-filled Kondo lattice 
model with attractive interactions. Typical truncation errors 
in the DMRG calculations are 5 x lO"'^ for U/J = -1.5. 
Errorbars are estimated from 1/L^ and 1/logL scalings. Gap 
energies, exchange constant J, and Coulomb interaction U are 
in units of t. 



First let us discuss the results for the case of strong ex- 
change coupling J = 10. Oi. The data for the excitation 
gaps are shown in Fig. ^. It is not difficult to identify 
the three phases as we scan U from to —2 J: the or- 
dinary spin liquid phase Ac > Ag > for U/J > —0.6, 
the spin liquid phase with identical spin and charge ex- 
citations Ac = As > for -0.8 > U/J > -1.4, and the 
gapless phase (SCDW) Ac = A^ = for U/J < -1.8. 
In the intermediate regime located between U/J = —0.8 
and —1.2 the difference of Ac and Ag vanishes in the 
bulk limit within the accuracy of the present calculation. 
Therefore this confirms the existence of the intermedi- 
ate phase where both spin and charge excitations have 
the same energy properties and are essentially of the two 
particle type as anticipated in the previous section. 

In addition the single quasiparticle excitation gap Aqp 
is shown in Fig. |[ The quasiparticles originate from 
doublons or holons and have therefore fermionic char- 
acter. The Aqp is obtained from the difference of the 
lowest energy in the Hilbert spaces with S^^^. — with 
N^ = L and S^^^ = ±1/2 with Nc = L±l. The Aqp is 
half of the Ac for U/J > —1.4. The reason is that the 
effective interaction between the quasiparticles is repul- 
sive in that region. On the other hand, for U/J < — 1.8 
the Aqp increases with growing ^U/J while Ac is zero 
(within our accuracy). Hence, we may interpret this as 
a switching from a repulsive to an attractive quasipar- 
ticle interaction, indicated by the minimum of the Aqp 
around U / J = —1.6. Since the character of the excitation 
changes here we expect that this switching coincides with 
the phase transition to the gapless phase. Obviously, in 
the thermodynamic limit all gaps should disappear at the 
transition point. Despite careful scaling analysis it is dif- 
ficult to determine the exact position of the transition 
from the present DMRG study because close to the tran- 
sition point the truncation error in the density matrix 
becomes large and the convergence of the RG iteration 
is rather slow. 

Next we turn to a weaker exchange coupling. The re- 
sults for J — 2.0t are shown in Fig. ||. In contrast to 
the strong coupling case we cannot find any indication of 
the intermediate regime with Ac = Ag > 0. This means 
that this regime is either absent or confined to a very 
small region close to the transition point which we lo- 
cate at around U/J « -2.0i. For U/J< -2.0t the two 
gaps essentially coincide and are practically zero. Cor- 
responding to this change the quasiparticle gap shows a 
minimum around the critical value of U. 

Similar calculations are carried out for J/t — 4,6,8 
and the phase diagram is obtained. In Fig. ^ we show 
the three phases by different dots which are numerically 
determined. The crosses in Fig. ^ represent estimated 
minimum points in the quasiparticle gap where we expect 
the phase transition to be located. There are technical 
limitations on the details of the phase diagram, which 
do not allow us to determine clearly the extension of the 
intermediate phase (Ac = Ag > 0). 



5 



12 



10'^ 



/ 6 

4 



Ac > As > 



□ □□□□□□ X 



Ac = As >0 



X □□□□□□-g 

Ac = As = 

□ □□□□□□□-C[] 



-10 

u 



-15 



-20 



FIG. 4. Ground state phase diagram of the half-filled 
Kondo lattice model with attractive interactions. The cross 
represents the point where the quasiparticle gap becomes min- 
imum. Exchange constant J and Coulomb interaction U are 
in units of t. 



sis shows that their correlation length is essentially the 
same. This means that not only the excitation gaps but 
also the velocities of the excitations are identical as we 
expect for the excitations of the particle-hole type. This 
clearly indicates that the picture of this phase as a kind 
of band insulator is appropriate. 
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A definite advantage of the DMRG scheme is that we 
can get very good approximations of the ground state 
wave functions for rather large systems. This allows the 
study of various correlation functions at least sufficiently 
far from the transition line of the phase diagram. Here 
we have analyzed the spin-spin and charge-charge cor- 
relation functions which show characteristic features of 
the phases. The correlation functions can be observed 
through the perturbing effect of the boundaries. This in- 
duces an oscillating disturbance into the wave function 
leading to a charge or spin density modulation analogous 
to the Friedel oscillations around an impurity. For exam- 
ple if Ac is finite then the density-density correlation has 
an exponential form and the charge density Friedel oscil- 
lations induced by an impurity potential shows the same 
exponential decay. The length scale is related to the ra- 
tio between the charge velocity and the gap, Wc/Ac. On 
the other hand, if Ac is zero, then a power law decay of 
the correlation function, similar to that of a Tomonaga- 
Luttinger liquid, is expected. This power law decay oc- 
curs naturally in the Friedel oscillations of both spin and 
charge. 

In Fig. ^ we show the charge and spin density Friedel 
oscillations, 5p{x) and <5(t(x), in the intermediate phase 
with J/t = 10 and U / J = 1.4. The charge density os- 
cillations are naturally induced by the open boundary 
conditions, while for generating the spin density oscilla- 
tions we have to apply a local magnetic field, Hiocai = 
2h{Sf — sf — Sl + sf^), coupling to the spins at both ends 
of the finite system. Obviously both the charge and spin 
density oscillations decay exponentially and our analy- 



FIG. 5. Charge density oscillations induced by the open 
boundary conditions, and spin density oscillations introduced 
by applying a local magnetic field. The system size is 30 sites, 
J — lot, and U = —lAt. The strength of the local magnetic 
field h is 0.2t. 

The charge and spin density Friedel oscillations in the 
gapless phase are shown in Fig. ^ and |^, respectively, 
by the solid line. In contrast to the exponential decay 
in the gapped phase, represented by the broken line, a 



power law decay is observed here, ^ {—ly/r 



Con- 



sidering the insets of Fig. ^ and ^ we find that the pow- 
ers 77c, s of decay are different for spin and charge (for 
the charge rjc « 1.15 and for the spin rjs « 2.4). This 
demonstrates the separation of spin and charge excita- 
tions as anticipated from the effective Hamiltionian in 
Eq.(7). In addition the staggering correlations both for 
spin and charge densities is obviously dominating. This 
is in agreement with the characterization of this phase 
as a spin and charge density wave with quasi long-range 
order, as the effective Hamiltonian in Eq.(7) suggests. 
From this result it is clear that among the two possi- 
bilities the Lieb-Schultz-Mattis theorem proposes for the 
phase —U:s>J the gapless state is the correct choice. 
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FIG. 6. Charge density oscillations of the half-filled Kondo 
lattice model with attractive interactions. The system size is 
30 sites and J — lOt. The solid line and the broken line 
correspond to U — —20t and — 14t, respectively. The inset 
shows log- log plot of the oscillations for U = —20t. 
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FIG. 7. Spin density oscillations of the half-filled Kondo 
lattice model with attractive interactions. The system size 
is 30 sites and J = lOt. The solid line and the broken line 
correspond toU — ~20t and — 14t, respectively. The strength 
of the local magnetic field h is 0.2t. The inset shows log-log 
plot of the oscillations for U = —20t. 



IV. CONCLUSION 

We have seen that the ID half-filled Kondo-Hubbard 
lattice with attractive electron-electron interaction has 
three phases, two gapped phases (spin liquid and band 



insultator-like) with short range correlations and one 
gapless phase with quasi long-range ordered spin and 
charge density waves (SCDW). Note that for positive U 
only the spin liquid exists. An indication of the differ- 
ence between positive and negative U can be found in the 
small- J limit. One may be tempted to lift the large spin 
degeneracy found for J = by perturbation theory, i.e. 
introducing the RKKY-interaction among the localized 
spins. In lowest order one separates the system into free 
conduction electrons and interacting localized spins. A 
simple calculation shows, however, that this perturbation 
concept fails for the ID system, because the effective in- 
teraction between the spins diverges for the wave vector 
q = 2kp = 77 for all U > 0. This indicates that there is no 
separation between these degrees of freedoms and numer- 
ical calculations suggest that for any finite J the ground 
state has spin liquid properties. In contrast, for negative 
U the perturbation converges for alLaiave vectors and one 
can derive a sensible RKKY-modelEj. Under this condi- 
tion the localized spins and the conduction electrons can 
separately exhibit gapless phases. Note, however, that 
the electron spin excitation has a gap for negative U and 
only the charge part is gapless with a tendency towards 
an enhancement of charge density wave correlations. The 
effective Hamiltonian (7) describes this kind of system in 
the extreme limit {\U\ ^ J). 

An interesting change occurs also for the nature of the 
excitations. If we scan from U = to [/ ^ J for fixed 
values of t and J then we find close to the transition 
{U « — 3J/2) a region where the lowest spin and charge 
excitations are due to doublon-holon excitation. There- 
fore, the two excitations have the same spectrum as a 
band insulator. We would like to point out, however, 
that the fermions, the holons and doublons, building the 
excitations are quasiparticles which are composed of the 
itinerant electronic and the localized spin degree of free- 
dom. For finite t they are real dressed quasiparticles 
involving complex correlation effects. Thus it is not an 
entirely trivial feature of this strongly correlated system 
to mimic the properties of a simple band insulator. On 
the other hand, in the spin liquid phase for \U\ <^ J the 
spin and charge excitation have different nature. The 
spin triplet excitation has excitonic character, i.e. it is a 
bound state of a holon-doublon excitation. 

Also in the gapless phase the spin and charge exci- 
tations are separated as is seen in the difference of the 
correlation functions. However, the origin of the spin- 
charge separation is different. The conduction electron 
and localized spins provide nearly independent degrees 
of freedom in close connection to the RKKY-model men- 
tioned above. 

The KHM with attractive interaction is an example 
of a strongly correlated electron system which possesses 
a number of phase of rather different character. These 
phases depend only on the relation of different coupling 
constants. Several of the mentioned features can be 
transfered to higher dimensional systems. The gapless 
ground state has very likely long-range order in both the 
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charge and the spin density wave correlation. However, 
it has to be noticed that the spin liquid phase for U > 
might have antiferromagnetic long-range order, if J is 
sufficiently small. This has not been proven so far. How- 
ever, this long-range ordered state has a charge gap (gap 
of a spin density wave state) and the spin excitations are 
due to both the electron and the localized spin. This is 
in contrast to the SCDW-state we discussed above. Thus 
in higher dimensions we may expect to see a richer phase 
diagram. 



2) |/, — /; S, S) : We consider the example |-|-, — ; TT) which 
leads to the matrix elements, 

mi = (+, -; TT |Weff|+, -; TT) = (-| + c)(^ + c') 
= (+,-;TT |Heff|-,+;TT) = ^(^ + c') 
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(A2) 

where the intermediate states consist of two triplets {E = 
—U + J/2) or of one singlet and one triplet {E = —U — 
J/2). 

3) |7, — /; S, —S): It is sufficient to consider the example 

l+,-;Ti), 



APPENDIX A: EFFECTIVE HAMILTONIAN IN 
THE LIMIT |Z7| > J > T 

The low-energy Hilbert subspace A in the limit of large 
attractive interaction consists of all real-space configu- 
rations of holons and doublons which appear in equal 
number at half-filling. All other states containing on- 
site singlets or triplets are higher in energy by a mul- 
tiple of order \U\. For convenience we introduced the 
notion of isospin I as an additional local SU(2) degree 
of freedom besides the localized spins. The doublon and 
the holon are mapped to isospin up and down, respec- 
tively (Eq.(6)). Thus each site has four basis states: 

\r,S^) = {|+,T>,|+,i),|-,T>,|-,i)}- For vanishing 
hopping t the whole subspace A is disconnected leading 
to complete degeneracy. We would now like to generate 
the efi'ective Hamiltonian which lifts this degeneracy and 
determines the low-energy physics. We use second order 
perturbation theory in the hopping where the higher en- 
ergy states outside of A appear as intermediate states. In 
this way the effective Hamiltonian will only have nearest 
neighbor coupling and be of the generic form, 

(Al) 

We consider now the coupling on a single bond, say be- 
tween site 1 and 2, where we represent the local states 
by |/f , Jf; 'S'l , /^l). There are three different cases to take 
into account: 

1) |/, /; S', its'): If the z-componcnt of the isospins on a 
bond are the same, then the hopping is ineffective and 
this state is not connected to any other via this bond. 



m2 = (+, -; Ti |Weff|+, -; Ti) = (-| + c)(-^ + c') 
= (+, -; Ti IWeffI-, +; Ti) = \{-^ + c') 



u-j u+j 



< 



(A3) 



and 



7713 



(+,-;Ti |Heff|+,-;iT) = (-^ + c)^ 
(+,-;Ti |Heff|-,+;iT) = 



(A4) 



2t^ 



3J 



> 0. 



Note that mi — m2 = ms. Now we can determine the 
coefficients in Eq.(7). Fixing a = 1 we find b = —1, 
c = -1/4, a' = b' and 



a' = 4(m2 - mi) < 
c' = — (mi -I- m2) > 0. 



(A5) 



Next, we rotate the isospin on every second site by tt 
around the z -axis, I^-y —> I^-^ = —I^-v and ^ = 
P by applying the phase factor exp(i7r/?). With this 
phase convention on the basis states b ^ —b and the 
effective Hamiltonian reaches its apparently SU(2) rota- 
tional symmetric form in both spin and isospin space as 
given in Eq.(7). 



8 



* On leave from Institute of Solid State Physics, University 
of Tokyo, Roppongi 7-22-1, Minato-ku, Tokyo 106, Japan. 

^ H. Tsunotsugu, Y. Hatsugai, K. Ueda and M. Sigrist, Phys. 
Rev. B46, 3175 (1992). 

^ T. Yanagisawa and K. Harigaya, Phys. Rev. B50, 9577 
(1994). 

3 A.M. Tsvelik, Phys. Rev. Lett. 72, 1048 (1994). 

* CC. Yu and S. White, Phys. Rev. Lett. 71, 3866 (1993). 

5 N. Shibata, T. Nishino, K. Ueda and C. Ishii, Phys. Rev. 

B53, R8828 (1996). 
® H. Tsunotsugu, M. Sigrist and K. Ueda, to be published in 

Rev. Mod. Phys. 

A. Auerbach, Interacting Electrons and Quantum Mag- 
netism, Springer, 1994. 
® H. Tsunetsugu, Phys. Rev. B55, 3042 (1997). 
E.H. Lieb, T.D. Schultz and D.C. Mattis, Ann. Phys. 16, 

407 (1961). 

I. Affleck and E.H. Lieb, Lett. Math. Phys. 12, 12 (1986). 
" S.R. White, Phys. Rev. Lett. 69, 2863 (1992). 

R. Egger and H. Schoeller, Phys. Rev. B54, 16337 (1996). 



9 



